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I. INTRODUCTION 

The degree to which the quark and gluon distribu¬ 
tions of the partons comprising the quark-gluon plasma 
(QGP) generated in relativistic heavy-ion collisions are 
momentum-space isotropic in the local rest frame (LRF) 
is currently an open question. There have been a num¬ 
ber of theoretical studies that have attempted to ad¬ 
dress this question using both perturbative QCD and the 
AdS/CFT framework (see Ref. [1] for a recent review). 
Ideally, however, one would like to have an experimental 
observable that could provide constraints on the degree 
of isotropy during the early stages of the QGP’s lifetime 
and perhaps, in addition, the subsequent approach to¬ 
wards isotropy. 

In principle, electromagnetic emissions are the ideal 
observable for studying the early-time dynamics of the 
QGP since they are weakly coupled to the plasma (a -C 
a s ). In addition, due to the fact that the QGP is ini¬ 
tially hot and then cools, high-energy (E > 2 GeV) pro¬ 
duction is dominated by early times when the system is 
in the QGP phase, while low-energy (E < 2 GeV) pro¬ 
duction receives significant contributions from late-time 
emissions when the system returns to the hadronic phase. 
This simple picture is complicated by the fact that there 
is a temperature distribution in the QGP, with the edges 
of the system being best described using hadronic degrees 
of freedom, however, since these regions are rather dilute 
and small in relative volume, the total radiation from 
this region is small compared to that produced from the 
central region. The two primary electromagnetic observ¬ 
ables studied in heavy-ion collisions are real photons and 
dileptons produced via decay of virtual photons. 

In this paper we focus on dilepton production from 


the deconfined phase of the QGP’s lifetime. The study 
of dilepton production from the QGP has a long his¬ 
tory, see e.g. Refs. [2-17]. For recent reviews, see also 
Refs. [18, 19]. Herein we focus on the effect of LRF 
momentum-space anisotropies on dilepton production. 
This work is an extension of previous studies performed 
in Refs. [20, 21] to include a realistic bulk evolution using 
the framework of anisotropic hydrodynamics [22-41]. For 
a recent review of the motivation for and methods used 
to obtain the anisotropic hydrodynamics equations and 
solve them numerically, we refer the reader to Ref. [42]. 

Herein, we make use of the anisotropic hydrodynamics 
equations obtained from the zeroth and first moments 
of the Boltzmann equation with the collisional kernel 
treated in the relaxation-time approximation and de¬ 
scribe the (3+l)-dimensional evolution of the QGP us¬ 
ing these equations. The resulting dynamical equations 
describe the full spatiotemporal evolution of the trans¬ 
verse temperature A and spheroidal momentum-space 
anisotropy parameter £ [43, 44]. The (3+1 [-dimensional 
framework allows both A, £, and the associated flow 
velocities to depend arbitrarily on the transverse coor¬ 
dinates, spatial rapidity, and longitudinal proper-time, 
however, herein we restrict ourselves to smooth Glauber- 
like initial conditions. 

The study presented herein is similar in spirit to prior 
studies of dilepton production using viscous hydrody¬ 
namics [15, 16]. In these works, however, the authors 
employed the standard viscous hydrodynamic lineariza¬ 
tion around an an isotropic thermal background. Our 
work goes beyond these studies by linearizing around 
anisotropic background and, as a result, we are able 
to better describe early-time dilepton production and 
dilepton production near the transverse and longitudinal 
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edges of the QGP. In addition, high-momentum dilepton 
production is treated in a more reliable manner since 
the anisotropic one-particle distribution function used 
to compute the dilepton rates is positive definite at all 
points in momentum space. We demonstrate that high- 
energy dilepton production is extremely sensitive to the 
assumed level of initial momentum-space anisotropy of 
the quark-gluon plasma. As a result, it may be possible 
to experimentally constrain the early-time momentum- 
space anisotropy of the quark-gluon plasma generated in 
relativistic heavy-ion collisions using high-energy dilep¬ 
ton yields. 

The structure of our paper is as follows. In Sec. II, 
we review the calculation of the leading-order dilepton 
production rate in an anisotropic QGP. In Sec. Ill we 
describe how one calculates the dilepton spectra includ¬ 
ing the effect of transverse and longitudinal expansion. 
In Sec. IV we present the setup for the anisotropic hy¬ 
drodynamics evolution, the resulting (3+l)-dimensional 
dynamical equations, the equation of state employed, and 
the initial conditions used. In Sec. V we present our final 
numerical results for the dilepton yields as a function of 
invariant mass, transverse momentum, and pair rapidity 
using fixed initial conditions and fixed final multiplicity. 
We present our conclusions and an outlook for the future 
in Sec. VI. We collect information about particle produc¬ 
tion from (3+l)-dimensional anisotropic hydrodynamics 
and compare this with Israel-Stewart viscous hydrody¬ 
namics in App. A. 


II. DILEPTON RATE IN ANISOTROPIC 
PLASMA 


We begin by reviewing the derivation of the dilepton 
emission rate for an anisotropic plasma starting from rel¬ 
ativistic kinetic theory. We follow the methodology pre¬ 
sented originally in Ref. [21]. The resulting formulas will 
be subsequently used in Sec. Ill to calculate the differ¬ 
ential dilepton spectra using anisotropic hydrodynamics 
framework. 

The dilepton emission rate is defined as the number 
of dilepton pairs produced per eight-dimensional phase- 
space volume 


dR l+l _ dN l+l 

d A P ~ d 4 Xd A P ’ ( ’ 

where = (t, x) and = (E , p) are the four-position 

and four-momentum, respectively. Based on relativistic 
kinetic theory , 1 at leading order in the electromagnetic 


coupling, <D(a 2 ), the dilepton emission rate follows from 


dR l+l 

d A P 


f d 3 pi d 3 p 2 , , 

J (2n) 3 (2tt) 3 •G (Pi) N ( p ' 2) 

X V qq a qq l £ ( 4 ) (P M -Pi -P2), 


( 2 ) 


where f q ^ is the phase-space distribution function of 
quarks (anti-quarks ), 2 v qq is the relative velocity between 
the quark and the anti-quark 


\/Pi • P2 rn 2 q 
2P p , 2P P2 


(3) 


and <7 l qq l is the total cross section for the leading-order 
quark-anti-quark annihilation process, q + q —> 7 * —»• 
l + +l~ 


i+r = / 2 m[\ / _ 4rn[\ 1/2 

99 3 M 2 { M 2 ) { M 2 ) 


(4) 


Henceforth, we will consider only high-energy dilepton 
pairs with invariant energies much greater than the lep¬ 
ton masses, M mi . Therefore, we will ignore lepton 
mass corrections appearing in Eq. (4) and simply take 
mi = 0 . 

Ultra-relativistic heavy-ion collisions are special in the 
sense that the matter created in such events is undergoing 
rapid expansion along the longitudinal (beam) direction. 
At the same time, the transverse expansion is initially rel¬ 
atively quite slow. One can show that this phenomenon 
inevitably leads to the presence of large momentum-space 
anisotropies in the phase-space distribution of the mat¬ 
ter. The simplest form for the distribution function that 
can be used to describe this situation is a generaliza¬ 
tion of an isotropic phase-space distribution which is 
squeezed or stretched along one direction in momentum 
space, defined by n, with a parameter —1 < £ < 00 , 
which describes the type and strength of the momentum- 
space anisotropy. In this case, the one-particle distri¬ 
bution function for the quarks and anti-quarks may be 
described at leading order by the following spheroidal 
“Romatschke-Strickland” form [43, 44] 

fg{q) (P. f. A ) = /?®(VP 2 + £(P-n) 2 , A), (5) 

where A is a transverse-momentum scale and £ is the 
anisotropy parameter introduced above. In the limit¬ 
ing case where £ = 0, Eq. (5) reduces to the standard 
isotropic distribution function. When £ = 0, A can be 
identified with the equilibrium temperature T of the sys¬ 
tem. Herein, we will take /* s ^ to be a Fernri-Dirac dis¬ 
tribution function f™?L(E,T) = [exp(£'/T) + 1] _1 . 

Using the Dirac delta function in Eq. (2), one can im¬ 
mediately perform the p 2 integration to obtain 


1 The same result can be obtained using standard finite tempera¬ 
ture field theory techniques. 


Prom now on we assume that fq = f q . 
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in (6) in the form 


dR l+l 5a 2 
d 4 P ~ 72tt 5 


d 3 pi 

E Pl E P2 


/ 9 (Pl:A,0/g(P2,A,0 


S(E-E P1 -E P2 ) = 2 (E- Pl )^j2 Wi - &i) > (8) 

VA , 


x 5(£-£ Pi -£ P2 ) 


P2=P-Pl 


( 6 ) 


To proceed, we parameterize the remaining three- 
momenta using spherical coordinates with the z-axis de¬ 
fined by the direction of anisotropy n, 


where 

X= (2p 1 Psin6 l p sinflpj 2 

— [2pi(E — PcosOp cosO Pl ) — M 2 } 2 . (9) 

The angles </>i are calculated as the two possible solutions 
to the equation 


Pi = pi (sin 9 P1 cos 4 > Pl , sin 9 Pl sin <j > pi , cos 9 Pl ), 

P = P(sin0pcos(/>p,sin0psin^>p,cos6 l p). (7) 

In this way we may rewrite the remaining delta function 


cos {<t>i - 4> Pl ) 


2p 1 (E — PcosOp cos 6 Pl ) — M 2 
2 piP sin 9 P sin 6 pi 


( 10 ) 


After these substitutions, we arrive at our final result for 
the dilepton emission rate 


dR l 


5a 2 


d 4 P 18tt 5 


J d(cos9 Pl ) J f q (p! yjl +£cos 2 9 P1 , 


c f q ^(P-pi) 2 +f(pi cos 9 Pl - P cos 9 P ) 2 , A^ 


( 11 ) 


with 

M 2 

ft± = 2{E — Pcos(0 P ± 9 Pl )) ' ^ ’ 

In order to evaluate the dilepton emission rate (11) it is 
necessary to perform the remaining two integrations nu¬ 
merically. In Fig. 1 we show the resulting dilepton emis¬ 
sion rate as a function of transverse momentum (left) 
and invariant mass (center), both scaled by A, and ra¬ 
pidity (right) for various values of anisotropy parame¬ 
ter £ S {—0.9,0,10,100} denoted by brown solid, red 
dashed, blue dotted and green dot-dashed lines, respec¬ 
tively. One can see that the production rate decreases 
(increases) due to increasing (decreasing) £. We note, 
however, that this is primarily due to the fact that in¬ 
creasing £ for fixed A results in a lower plasma density. 
In order, to properly assess the impact of anisotropies on 
the production, one has to fold these rates together with 
a realistic model of the full spatiotemporal evolution of 
both £ and A. 

III. DILEPTON SPECTRA 

Our final goal is to study the impact of space-time 
dependent anisotropies in the system on the dilepton dif¬ 
ferential spectra. In this way we hope to probe the early 
stages of the quark-gluon plasma, where the anisotropies 
are expected to be the largest. In order to do this one 
must include in Eq. (11) the space-time dependence of 


A and £ using some hydrodynamic model and then inte¬ 
grate over the entire space-time volume (which contains 
the quark-gluon plasma phase) and the appropriate mo¬ 
menta/invariant mass cuts for the dilepton pairs. For 
this purpose, we parametrize the pair four-momentum in 
the standard way, 

= (m_L coshy, p± cos 4> p ,p± sin</> p , m± sinhy), (13) 

where to_l = \JM 2 + p ■} defines the transverse mass and 
y = 1/2 ln[(S + P ||)/(E — py)] is the momentum-space 
rapidity. Above, we used p ±, p\\, and <p p to denote trans¬ 
verse momentum, longitudinal momentum, and momen¬ 
tum azimuthal angle, respectively. 

One can also use the usual Milne hyperbolic 
parametrization of space-time which is convenient for de¬ 
scribing heavy-ion collisions within the relativistic hydro¬ 
dynamics framework 

= (TCOshc,xj_, r sinh<r). (14) 

In Eq. (14), we used r = \ft 2 — z 1 and c = tanh _1 (z/t) to 
denote the longitudinal proper time and the space-time 
rapidity, respectively. With these parameterizations, the 
differential measures for four-momentum and space-time 
are d 4 P = MdM dyp±dp±dcj) p and d 4 X = rdr d<; d 2 x±, 
respectively. This allows us to calculate the invariant 
mass and transverse momentum differential spectra using 


















4 




FIG. 1. The dilepton emission rate as a function of transverse momentum (left), invariant mass (center) and rapidity (right). 
For the transverse momentum dependence (left) we fixed M/A= 3 and y= 0, for the invariant mass dependence (center) we 
fixed p_ l/A= 3 and y= 0, for rapidity dependence we fixed p ±/A = M/ A= 3. 


dN l+l 
M dM dy 

dN l+l ~ 

p±dp±dy 


/ p±dp± 

/pmin 

»M max 

/ MdM 
I M min J 


d(j) p fdrX 


dR l+r 

d 4 P 


d(j) p d A X 


dR l+l ~ 

d A P 


(15a) 

(15b) 


respectively, where the integration ranges p± ln , p™ ax and 
M mm , M max will be specified later according to the ap¬ 
propriate physical/experimental cuts. The integration 
over the space-time volume is performed only in the de- 
confined quark-gluon plasma stage. In practice, we only 
include contributions from regions that have an effective 
temperature that is higher than a critical temperature, 
i.e. T = 7?. 1 / 4 (^)A > T c with 7 Z(£) defined in Eq. (45). 
In all results shown herein, we assume T c = 175 MeV. 


We will assume that when the system reaches T c , all 
medium emission stops. We do not take into account the 
emission from the mixed/hadronic phase at late times 
since the kinematic regime we study (high M and p±) is 
dominated by early-time high-energy dilepton emission. 
Due to the large uncertainty connected with the correct 
value of critical temperature existing in the literature we 
also checked that the results obtained here are almost 
completely independent of the choice of T r in the range 
150 - 200 MeV. 

Equations (15a) and (15b) are evaluated in the cen¬ 
ter of mass of the colliding nuclei (LAB) frame while the 
dilepton emission rate is calculated in the local rest frame 
(LRF) of the emitting region. Therefore, before evaluat¬ 
ing Eqs. (15) we have to boost the LAB frame momentum 
to the LRF of the fluid cell using p= A^ v p v , where 
the Lorentz boost tensor 


( 7 

-7W 

-7 Vy 

-7Az \ 

-7W 

l + (7-l)$ 

(7-1)^ 

(7-1)^ 

-7 Vy 

(7-1)^ 

l + (7-l)S 

(7-1)^ 

\ ~JVz 

(7-1)^ 

(7-1)^ 

l + (7-l)ui / 


(16) 


depends on the four-velocity of the fluid element 
u M ( x M ) = 7(1 ,v x ,v y ,v z ), where 7 = 1/y/l — v 2 and 

v = \Jv1 + Vy + v\ . One can easily check that, as ex¬ 
pected, Uj)rf = ul/ = (1,0,0,0). Making use of 
Eq. (11) in Eqs. (15), we obtain the dilepton spectra in¬ 
cluding the effect of a space-time-dependent momentum 
anisotropy. 


IV. HYDRODYNAMIC EVOLUTION 

As mentioned above, in order to make predictions for 
the differential dilepton spectra expected to be produced 
from the QGP phase, one must integrate over the full 


space-time history of the QGP. For this purpose, we use 
anisotropic hydrodynamics. Anisotropic hydrodynamics 
reduces to second-order viscous hydrodynamics in the 
limit of small anisotropy [41], but reproduces the dy¬ 
namics of the QGP more reliably when there are large 
momentum-space anisotropies. 

A. (3+l)-dimensional anisotropic hydrodynamics 

In this paper, we assume that the system created dur¬ 
ing the collision of the heavy ions evolves through a non¬ 
equilibrium state and that the quark and anti-quark one- 
particle distribution functions are well approximated by 
Eq. (5) both at early times and late times. At the same 
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time, we assume that, although the system is highly 
anisotropic it may still be, to good approximation, de¬ 
scribed using hydrodynamic-like degrees of freedom, such 
as energy density and pressures. 3 In this way, the de¬ 
tailed microscopic description of the system can be re¬ 
placed by an effective description which realizes simple 
physical laws, such as conservation of energy and mo¬ 
mentum. In the following, we will present the frame¬ 
work of leading-order anisotropic hydrodynamics [22- 
26, 28, 29] which is designed to describe a potentially 
highly-anisotropic plasma by assuming that its distribu¬ 
tion function is, to good approximation, expressible in 
the form given by Eq. (5). 

At leading order, one can derive the equations of mo¬ 
tion of the anisotropic system starting from kinetic the¬ 
ory assuming that the distribution function of the sys¬ 
tem is known. This can be done by taking moments of 
the Boltzmann kinetic equation with the collision term 
treated in the relaxation-time approximation (RTA) 

P fl d ll f = P ^(f iso -f), (17) 

Teq 

where r eq is the microscopic relaxation time which can 
depend on position and time. Taking the first moment of 
the Boltzmann equation results in the energy-momentum 
conservation equation 


= 0. (18) 

Taking the zeroth moment of the Boltzmann equation 
results in the particle production equation 


Equations (18) and (19) provide a set of five indepen¬ 
dent partial differential equations 


D u s — — (e + P±) 9 U + (Pj_ — Py) u^DzZ 11 , 
D Z P\\ = {P± — P\\) 9 Z + {e + P±) z v D u u u , 


D u u± 


u _L U_L ■ Vj.Pl 

£ + P± u]_ 


+D U P± + (Pj l — P\\)u v D z z 1 ' 



1 

y-l(£ + P±) 


(j-ix^y 


U y 3 X ) P _L , 


and 


( 22 ) 

(23) 


(24) 

(25) 


D u i 3D U A 

2(1 + 0 A 


9 U ~\ - 

7"eq 


i-p 3/4 (0v / TT£ 


(26) 


respectively, for five parameters: the four-velocity vf 1 , the 
transverse temperature A, and the anisotropy parameter 
£. 4 * In the above equations, we use a + subscript to indi¬ 
cate two-dimensional vectors in the transverse plane, e.g. 
uj_ = ( u x ,u y ) and V_l = ( d x ,d y ). We have also intro¬ 
duced a compact notation for the convective derivative 
D u = u^d^, the longitudinal derivative D z = z^d^, and 
the expansion scalars d u = and 9 Z = d^z 11 . 

In the most general case, where the matter expands 
in the longitudinal and transverse directions without 
any symmetry constraints, one can use the following 
parametrization of the LAB frame four-velocity of the 
fluid and the space-like four-vector z M 


dpN 1 * = 


'T’eq 


(19) 


At leading order, the energy-momentum tensor has the 
form typical for a spheroidally anisotropic system 


= (e + P ± ) «V - P ± - (P ± - P[\)z lt z v , (20) 

and the particle flux is defined in the standard manner 

A+ =n eq ?+. (21) 

In Eqs. (20) and (21) e, n, P \\, and P± stand for energy 
density, particle density, longitudinal pressure, and trans¬ 
verse pressure, respectively. The four-vector z M is orthog¬ 
onal to u M and in the LRF points in the longitudinal di¬ 
rection (identified with the direction of the anisotropy in 
the system, n) [28]. 


3 This assumption has been tested elsewhere by comparing the 
predictions of anisotropic hydrodynamics to exact solutions of 
the Boltzmann equation in a variety of special cases [33, 34, 
40, 45-49]. These studies found that anisotropic hydrodynamics 
provides the most accurate description of both the early and late 
time behavior of QGP dynamics. 


= (uocosh'i?, uj_,uosinh$), (27) 

z* = (sinht?,0,coshd), (28) 

where we introduced the longitudinal rapidity of the fluid 
cell '(). Using the four-velocity normalization condition, 


td'iq, = 1, one has 


u ± = \J u l + u l ■ ( 29 ) 

With the parametrizations (27) and (28), one may calcu¬ 
late the following quantities appearing in Eqs. (22)-(26), 

D u = u_l • V_l + uqLi , (30) 

9u = V_l ■ u_L +LiUo + u 0 L 2 'd, (31) 

D z =L 2 , (32) 

9 Z = L rtf, (33) 

u v D z z v = uoL-2'9 , (34) 

z u D u u v = -u 0 (u_l • Vj_ +u 0 Li) i?, (35) 


4 Note that the four-velocity satisfies = 1 and hence it con- 

tains only three independent degrees of freedom. 
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where the two linear differential operators, L\ and L 2 , 
are given by 


Li = cosh(c — d)d T — sinh(c — 1 ?)—, (36) 

T 

—1/2 = sinh(c — $)3 r — cosh(c — i?) —. (37) 

T 


We also use the relation between the relaxation time r eq 
and the shear viscosity to entropy density ratio fj = i)/s 
[ 22], 5 


5?y 

Teq= 2 T' 


(38) 


B. Anisotropic equation of state 


Herein, we consider a system that consists of massless 
particles described by the anisotropic distribution func¬ 
tion (5). Using standard kinetic theory definitions 


N^ = J d 3 PpVf, 

(39) 

T ^ = J d 3 Pp»p v f 1 

(40) 


where d 3 P = d 3 p/ [(27r) 3 p°], and the tensor decomposi¬ 
tions specified in Eqs. (20) and (21), one can calculate 
the thermodynamic properties of the system 


IK C\ n ' ls ° (A) 
n(A.t) = . , 

v Vi + C 

(41) 

£ (A, 0 = 72.(0 £ iso(A), 

(42) 

-Pj_(A, £) = 7£_l(£) -Pi SC) (A), 

(43) 

P\\ (A, 0 = 7£|| (0 -Piso(A), 

(44) 


where ?ti so , £i so , and -Piso are the isotropic particle density, 
energy density, and pressure, respectively, and 


S 2 

1 tan 1 

i + £ 




£ +1 


^ii(0 = | 

K + W-i' 
e + i 



(45) 

(46) 

(47) 


Herein, we assume the simple case of a conformal fluid, 
i.e. £j so = 3Pj so . As a result, Eqs. (41)-(44) describe 
the equation of state of an anisotropic system of classical 
massless particles with vanishing chemical potential. 


5 We note that the factor of 2 in the denominator of Eq. (38) 
is needed if one uses the spheroidal form (5) together with the 
zeroth and first moments of the Boltzmann equation. 


C. Initial conditions 


In order to solve the set of partial differential equa¬ 
tions (22)-(26) in general (non-boost-invariant (3+1)- 
dimensional evolution), one has to make a reasonable 
assumption about the initial conditions at the initial lon¬ 
gitudinal proper-time for the hydrodynamic evolution, 
r = To, i.e. one has to define five three-dimensional pro¬ 
files: A(r 0 ,xj_,<r), £(t 0 ,x_ l ,c), Mx(t 0 ,x_l,<j), %(t 0 , xj_, c), 
and $(t 0 ,xj_,c). 

During a heavy-ion collision, due to inelastic interac¬ 
tions the participating nucleons deposit some energy in 
the space-time volume of the fireball. In this work, we 
assume that the distribution of deposited energy is well 
described by the optical Glauber model. 6 Herein, we as¬ 
sume that the initial energy density is proportional to 
the scaled initial density of the sources. Therefore, the 
transverse momentum scale is given by 


A(t 0 ,xj_,?) =£ iso 1 


’ P{b,x_ L,c) \ 

e ° p(o> 0 ,0) ) ’ 


(48) 


where the proportionality constant £0 is chosen in such a 
way as to reproduce the total number of charged particles 
measured in the experiment, and £^ denotes the inverse 
£iso(A) function. 7 

The density of sources is constructed using the follow¬ 
ing mixed model 


P(b,x_ L ,c) 


(! - ^(PwN^+iJ+PwNtMi)) 


+ 2n p-Qc{b, xj_) 


/(? - ?s(&,xjJ), (49) 


where is the density of wounded nucleons from the 
left/right-moving nuclei and psc is the density of binary 
collisions, both of which are obtained using the optical 
limit of the Glauber model 


Pwn(^ X i) — ^ ^ X -L+ 2 ^ 


1 —e 




PBc{b,X±) = <JinT 



T 



,(50) 

(51) 


The longitudinal profile is taken to be 

(c - Ac) 2 


/(c) = exp 


2 <U 2 


©del-Ac) 


(52) 


For the LHC case studied here, we use n = 0.145 for 
the mixing factor and an inelastic cross-section of cr,; Tt = 


6 Although it is quite interesting, for this first study we do not take 
into account initial fluctuations in the position of the nucleons 
or nucleonic substructure. We postpone the Monte-Carlo event- 
by-event analysis to a future work. 

7 In principle, one could use the full expression for the energy 
density given by Eq. (42) in Eq. (48), however, since we use 
an initial anisotropy profile that is homogeneous in space, this 
would merely result in the overall multiplicative factor which can 
be absorbed by rescaling £ q . 
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62 mb. We also restrict ourselves to the minimum-bias 
studies in which case b = |b| = 9.5 fm. The parameters 
of the longitudinal profile (52) were fitted to reproduce 
the pseudorapidity distribution of charged particles with 
the results being Ac = 2.5 and <r ? = 1.4. The shift in 
rapidity is calculated according to the formula [50] 


_ _ 1 l„ %N + %N + Pwn) /con 

SS = x in — ; - rq: = r > l 0 ' 3 ; 

Pwn + Pwn v P\P\vn Pwn) 


where all functions are understood to be evaluated at a 
particular value of b and x_l. The participant velocity is 
defined as vp = s/{\/s/ 2) 2 — (tojv/2) 2 /(v / s/2) and mjy 
is the nucleon mass. In Eqs. (50)-(51) we have made use 
of the thickness function 

T{x±) = j dz pws(xi,z), (54) 


where the nuclear density is given by the Woods-Saxon 
profile 


Pws ( X _L, z) = Po 


1 + exp 


^V X J. 2 + z 2 ~ RJ 


(55) 


For Pb-Pb collisions, we use po = 0.17 fm -3 for the nu¬ 
clear saturation density, R = 6.48 fm for the nuclear 
radius, and a = 0.535 fm for the surface diffuseness of 
the nucleus. 

In the calculations presented herein, we assumed that 
the produced matter has initially no transverse flow, i.e. 
u x (r 0 ,xj_,c) = u x (t 0 , xj_, c) = 0, while the initial lon¬ 
gitudinal flow is of Bjorken form $(to,xj_,c) = C- For 
simplicity, the initial anisotropy parameter is assumed to 
be homogeneous, £(t 0 ,xj_,c) = Co - 8 


four-dimensional interpolating function from the table, 
we could assume that it is valid at continuous values of 
these four variables. We then evaluated the remaining in¬ 
tegrations over space-time, transverse momentum angle, 
and transverse momentum or invariant mass appearing 
in Eqs. (15) using Monte Carlo integration. For the inte¬ 
gration over the transverse momentum we have specified 
the default cuts as follows: p“ ln = 1 GeV and p“ ax = 20 
GeV, while for the invariant mass integration we used 
M min = 1 GeV and M max = 20 GeV. 


A. Dilepton production with fixed initial conditions 

We begin by presenting the minimum-bias dilepton 
spectra at midrapidity, y = 0, calculated assuming fixed 
initial conditions. The initial central temperature was 
taken to be T 0 = 567 MeV, at a starting time of r 0 = 0.3 
fm/c. For this case, we fixed the initial anisotropy pa¬ 
rameter to Co = 0, which means that the system is ini¬ 
tially isotropic in momentum space. In Fig. 2 we plot 
the resulting invariant mass spectra (left) and transverse 
momentum spectra (right) of dilepton pairs for various 
values of fj. We can see that the spectra flatten with in¬ 
creasing fj and the normalization increases with increas¬ 
ing fj. The latter implies that the total final multiplicity 
changes with changing fj. This effect is more visible in the 
p _l spectra. We note here that the fixed-initial-condition 
behavior described above is monotonic in fj. The increase 
in multiplicity with increasing fj is related to dissipative 
particle production in the QGP, which would also be re¬ 
flected in increased final particle multiplicity across all 
particle types. 


V. RESULTS 

In this section, we present our model predictions for 
the minimum-bias e + e _ yields resulting from Pb-Pb col¬ 
lisions at LHC with yfs = 2.76 TeV beam energy. Before 
presenting our results, we first explain the setup and pa¬ 
rameters chosen for our calculations. 

Since the differential clilepton rate dR l 1 /d A P given in 
Eq. (11) is independent of the assumed space-time model, 
we first evaluate it numerically using double-exponential 
integration on a uniformly-spaced 4-dimensional grid in 
M/A, p±/A, y, and log 10 (£ + 1) such that M/A,p±/A € 
{0.1,40}, y G {-6,6} and log 10 (£ + 1) £ {-1,3}. 9 The 
spacing was chosen in such a way that, after building a 


8 On general grounds, one can expect that the level of momentum- 
space anisotropy is larger in regions that have a lower effective 
temperature. As a result, our assumption of a constant £o is a 
conservative one. 

9 Note that the dilepton rate is an even function of y. therefore, in 
practice, we may restrict ourself to positive values of y only. 


B. Dilepton production with fixed final multiplicity 

The case presented in Section V A is unphysical since 
the average final particle multiplicity in a given centrality 
range is held fixed when presenting experimental results 
for the dilepton spectra. We presented the prior case 
only to establish that, for fixed initial temperature, the 
behavior seen in the final dilepton spectra is monotonic 
when fj is increased. In this section, we present the same 
analysis, however, now, for each presented case, the ini¬ 
tial central temperature is rescaled in such a way as to 
keep the final multiplicity of particles at freeze-out fixed. 

Fig. 3 presents the invariant mass spectra (left) and 
transverse momentum spectra (right) of dilepton pairs 
at midrapidity, y = 0. In the left panel, we can clearly 
see the effect of the rescaling of the initial temperature 
for cases with various fj, i.e. the spectra does not change 
significantly as long as the final multiplicity of particles 
at freeze-out is fixed. However, importantly, we observe 
that the spectra do not necessarily have a monotonic de¬ 
pendence as fj is increased (see e.g. the left panel of 
Fig. 2). 
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FIG. 2. (Color online) The invariant mass spectra (left) and transverse momentum spectra (right) of dilepton pairs at 
midrapidity, y = 0, for various values of shear viscosity to entropy density ratio 4nfj £ {0,0.1,1,2, 3, 6,10}. For all cases the 
initial temperature is fixed to To = 567 MeV and the system is initially isotropic in momentum space, £o = 0. 




FIG. 3. (Color online) Same as Fig. 2, except here, instead of fixing the initial temperature, we keep the final particle 
multiplicity fixed. 


The non-monotonic behavior primarily due to the fact 
that particle production within anisotropic hydrodynam¬ 
ics is not a monotonic function of fj , as it is for standard 
viscous hydrodynamics [28, 33, 34]. Instead, one observes 
a maximum in particle production at a certain value of fj 
which depends on the assumed initial temperature. The 
fact that there must be a maximum can be anticipated 
by the fact that particle production should vanish in the 
both the ideal and free streaming limits. 10 As a result, 
when fixing the initial temperature to guarantee fixed 
final multiplicity, the required temperature may not be 


10 The dependence of particle production on the assumed value of 
fj is discussed in more detail in Appendix A. In that appendix, 
we compare particle production as a function of fj using both 
anisotropic and viscous hydrodynamics. 


monotonically decreasing as ij is increased. 


Note that, although the final multiplicity of particles 
created at freeze-out is fixed, the number of dileptons 
which are produced in the QGP volume varies with fj. As 
a result, we observe a small but noticeable non-monotonic 
change in the dilepton invariant mass spectra. Similar ar¬ 
guments also apply to the transverse momentum spectra 
shown in the right panel of Fig. 3. For the p_i_-spectra 
the effect is smaller. Based on our final results shown in 
Fig. 3 one can see that for Anfj £ (1, 3), which spans the 
range of fj extracted from the flow experimental data, the 
impact of shear viscosity in the system on the dilepton 
spectra is quite small. 
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FIG. 4. (Color online) The invariant mass spectra (left) and transverse momentum spectra (right) of dilepton pairs at 
midrapidity, y = 0, for various initial anisotropy conditions, and 4nfj = 1. The results with £o = —0.9,0,10, and 100 are 
denoted by brown solid, red dashed, blue dotted and green dot-dashed lines, respectively. 




[GeV] 


FIG. 5. (Color online) Same as Fig. 4 however here we take 4-7 t?7 = 3. 


C. Effect of initial anisotropy 

We now turn to the analysis of the impact of the ini¬ 
tial anisotropy in the system, £g, on the dilepton spectra. 
In Figs. 4 and 5 we present invariant mass spectra (left 
panels) and transverse momentum spectra (right panels) 
for dilepton pairs at midrapidity, y = 0 , for various initial 
anisotropy conditions. The results with £ 0 = —0.9,0,10, 
and 100 are denoted by brown solid, red dashed, blue 
dotted and green dot-dashed lines, respectively. The val¬ 
ues of £o < 0 (£o > 0 ) correspond to prolate (oblate) 
initial momentum distribution functions. In Fig. 4 we 
keep the viscosity fixed to fj = 1/A.it while in Fig. 5 we 
set fj = 3/4-7T. In each case the initial energy density at 
the center, eg, is rescaled to keep the final multiplicity 
of particles at freeze-out fixed. The values of £q used in 
each case are listed in Table I. From Figs. 4 and 5, one 
can clearly see that the transverse momentum spectra are 


47t?7/s^\^^ 

-0.9 0 10 100 

0.1 

- 72.11 

1 

12.98 64.53 235.86 714.31 

3 

13.69 60.44 215.61 660.74 


TABLE I. Values of the initial central energy density, 
£q [GeV/fm 3 ], used in all the figures of this Section except 
for Figs. 2 and 3. 

quite sensitive to the initial momentum anisotropy in the 
system. For an initially oblate configuration, they are be¬ 
coming flatter. The opposite behavior is observed for an 
initially prolate configuration. This effect is particularly 
significant for large values of p±. This opens possibility 
to measure initial anisotropy of the plasma by looking at 
large p± dilepton pairs at LHC. The behavior of the in- 
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FIG. 6. (Color online) The rapidity dependence of the dilep¬ 
ton modification factor &(fj) for 47T7) £ {0.1,1,3} denoted by 
red dashed, blue dotted and green dashed-dotted lines, re¬ 
spectively. The initial anisotropy is £o = 0 in this case. In 
this case we also use default cuts: p± tn = 1 GeV, p™ ax = 20 
GeV, M min = 1 GeV and M max = 20 GeV. 

variant mass spectra, on the other hand, is more difficult 
to understand since, in this case, both oblate and prolate 
initial conditions lead to a flattening of the spectra. 

D. Production at forward rapidities 



y 


We close this Section by presenting an analysis of dilep¬ 
ton production at the forward rapidities following the 
preliminary study made in Ref. [51]. In Fig. 6 we present 
the dilepton modification factor 

(56) 

for 47 t? 7 € {0.1,1,3} and fixed value of £o = 0. In this 
figure, one sees that increasing the value of fj results in a 
suppression of particle production. The emission is more 
suppressed when we go to more forward rapidities, up to 
40% in the most extreme case. In Fig. 7, we present a 
complementary study of the dilepton modification factor 
(analogous to (56)) 


*(£<>) 


dN* + e-(Z q) ^ j V e+e -(6) = 0 ) ^ _ (57) 


for £ 0 = —0.9, 0,10 and 100 (the notation is the same as 
in Section V C) and for Anfj = 1 (top panel) and 47 t? 7 = 3 
(bottom panel). Similarly to Fig. 6, we observe a sup¬ 
pression of the dilcpton production at forward rapidities, 
which increases with increasing initial anisotropy param¬ 
eter £ 0 and the viscosity in the system. Moreover, we 
observe the opposite effect when the distribution is ini¬ 
tially prolate. In this case, we find dilepton enhance¬ 


FIG. 7. (Color online) The rapidity dependence of the dilep¬ 
ton modification factor 4>(£o) for £o = —0.9,0,10 and 100 
(notation is the same as in Section V C) and for 47 t? 7 = 1 (top 
panel) and 4-7 t? 7 = 3 (bottom panel). The p± and M cuts are 
the same as in Fig. 6. 


ment at forward rapidities together with moderate sup¬ 
pression in midrapidity. These effects provide the pos¬ 
sibility to probe the initial degree of thermalization of 
the system by looking at forward rapidity emission of 
dilepton pairs. Finally, in Fig. 8, we present the dilepton 
enhancement factor $ with different cuts corresponding 
to 3 GeV < p ± < 20 GeV and 3 GeV < M < 20 GeV. As 
can be seen from this figure, high-energy dilepton emis¬ 
sions are more sensitive to the level of momentum-space 
anisotropic in the quark-gluon plasma. Of course, since 
statistics are more limited, high-energy dilepton spectra 
are usually more difficult to measure accurately. 


VI. CONCLUSIONS 

In this paper we computed the dilepton invariant 
mass and transverse momentum spectra produced from 
the quark-gluon plasma. To accomplish this, we used 
the leading-order (3+l)-dimensional anisotropic hydro- 
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FIG. 8. (Color online) Same as the bottom panel of Fig. 7, 
except high-energy cuts of 3 GeV < p± < 20 GeV and 3 GeV 
< M < 20 GeV. 


dynamics equations obtained from the zeroth and first 
moments of the Boltzmann equation and assumed a con¬ 
formal (ideal) equation of state. The anisotropic hydro¬ 
dynamics equations solved allow for both azimuthal spa¬ 
tial anisotropy and a realistic rapidity profile. In this pa¬ 
per we considered a fixed (min-bias) impact parameter. 
We found that, when adjusting the initial temperature in 
order to enforce fixed final particle multiplicity, both the 
dilepton invariant mass and transverse momentum spec¬ 
tra show only a weak dependence on the assumed value 
of 77 /s. 

A similar conclusion was found in an earlier works 
that used a much more primitive model of the dy¬ 
namics [21, 51]. With the inclusion of the full (3+1)- 
dimensional dynamics using anisotropic hydrodynamics, 
we are now more confident that the dilepton spectra only 
have a weak dependence on the assumed value of 77 /s. 
That being said, in these previous works the possibil¬ 
ity of a finite initial momentum-space anisotropy £0 was 
not considered. In this work we found that the high- 
mass and high-transverse-momentum dilepton spectra 
are quite sensitive to the initial level of momentum-space 
anisotropy. Additionally, we demonstrated that the ra¬ 
pidity dependence of dilepton production is also sensitive 
to the initial level of momentum-space anisotropy. These 
observations offer some hope that one might be able 
to experimentally determine information about early- 
time momentum-space anisotropies generated in heavy- 
ion collisions using dilepton production. 

In this work we made a few simplifying assumptions 
that will be improved in future works. The first of these 
is that we only study min-bias collisions. The magni¬ 
tude of the effects seen here could depend on centrality 
in a non-trivial way since in central collisions the plasma 
lifetime is significantly longer but the level of momentum- 
space anisotropy developed dynamically in the center of 
the fireball will be reduced. We plan to make a sys¬ 


tematic study of the centrality dependence of our re¬ 
sults in a forthcoming paper. Another crucial assump¬ 
tion was that we used only the leading order (Born) rate 
for dilepton production. It is possible that inclusion of 
the next-to-leading order rate could significantly mod¬ 
ify our conclusions. Unfortunately, to the best of our 
knowledge such a calculation only exists for an isotropic 
quark-gluon plasma [13]. It would be very interesting to 
see if these calculations could be extended to the case of 
an anisotropic quark-gluon plasma. 

Looking forward, one should also consider polarized 
dilepton emission as suggested in Ref. [52]. The polar¬ 
ization asymmetry could be quite sensitive to early-time 
momentum-space anisotropies and possibly also to the 
assumed value of r)/s. Finally, we mention that another 
ideal observable that should be studied further is the 
emission of real photons. This has been studied using 
viscous hydrodynamics in Refs. [53-55] and using simple 
models of anisotropy evolution in the plasma [56-62], It 
is necessary to extend these studies to include the (3+1)- 
dimensional evolution of the QGP using anisotropic hy¬ 
drodynamics in order to draw more firm conclusions 
about the effect of momentum-space anisotropies on pho¬ 
ton production. We also mention that, like dileptons, 
a difficult, but necessary, step will be to extend the 
NLO calculation of photon production first obtained in 
Refs. [11, 12, 14] to an anisotropic quark-gluon plasma. 
The difficulty in this calculation stems from the presence 
of color plasma instabilities that render the NLO rate 
formally infinite. In practice, these infinities will be reg¬ 
ulated due to the eventual saturation of unstable mode 
growth, but how to implement this in practice is an open 
question. 
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Appendix A: Particle production in viscous and 
anisotropic hydrodynamics 

One can show that particle production within 
anisotropic hydrodynamics is not a monotonically in¬ 
creasing function of shear viscosity to entropy density, 
77 /s. This behavior is in agreement with exact solutions of 
RTA Boltzmann equation [28, 33, 34]. The fact that there 
must be a maximum in particle production as a function 
of 77 /s can be anticipated by the fact that particle produc¬ 
tion should vanish in both the ideal and free streaming 
limits. The behavior found using anisotropic hydrody¬ 
namics is qualitatively different than all known standard 
second-order viscous hydrodynamics approaches, which 
predict that particle production increases monotonically 
as 77 /s increases. The non-monotonicity of particle pro- 
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duction becomes particularly important when enforcing 
fixed final multiplicity of particles, since this is typically 
accomplished by rescaling the initial central temperature 
while holding other parameters fixed. Such a temper¬ 
ature rescaling can affect dilepton yields, since there is 
strong sensitivity of the dilepton spectra to the temper¬ 
ature of the emitting source. 

In order to extract the freeze-out hypersurface, we pa¬ 
rameterize space-time in the following way 

t = (r 0 + d( C, </>, 9) sin 8 sin Q cosh(d(£, <j>, 8) cos 8 ), 
x = d(C, </>, 8) sin 8 cos ( cos <j ), 
y = d(£, (f>, 8) sin 8 cos ( sin (j ), 

2 = (to + d(C, 0, 9) sin 8 sin £) sinh(d(£, </>, 8) cos 8 ). (Al) 

This parametrization leads to simple formulas for the 
space-time rapidity longitudinal proper time r, and 
the transverse distance r, 

? = d(C, </>, 8) cos 8 , 

t — To T d(C, 4>, 9) sin 8 sin (, 

r = d(C, 8) cost). (A2) 

The three angles £, <j>, and 8 are restricted to the ranges 

0 <C< tt/2, 

0 <4>< 2n, 

0 < 8 < 7T. (A3) 

The quantity d((,(j),6) describes the distance between 
a point on the freeze-out hypersurface and the coordi¬ 
nate system’s origin (r = r 0 ,x = 0 ,y = 0,<r = 0). The 
parametrization (Al) works quite well for all smooth ini¬ 
tial conditions where the distance d is a function of £, 
(j), and 8. Using the parametrization (Al), one can inte¬ 
grate the particle number on the freeze-out hypersurface 
specified by constant effective temperature Xpo i n the 
following way 

N = f dSX^(?FO^- 1/4 (e(C,^0)),e(C^^)) , 

( A4 ) 

where the form of dT,^ may be obtained with the help of 
the formula known from differential geometry 

d^=e, a0 , — — —d(dW8. (AS) 


The tensor e^ap-y is the four-index antisymmetric Levi- 
Civita tensor with £0123 = 1. The quantity defines 
a four-vector that is perpendicular to the hypersurface 
at point x^. Its norm is equal to the volume of the hy¬ 
persurface element. The variables £, <j>, and 8 introduce 
a coordinate system in Minkowski space parameterizing 
the positions of points on the freeze-out hypersurface. 
Their ordering is chosen in such a way that points 
in the direction of decreasing temperature. 

In Fig. 9, we plot the particle production measure, 
Naniso/Aideal-1 as a function of fj, where N aniso (ideal) de¬ 
notes the density of gluons (41) integrated on the isother¬ 
mal hypersurface, i.e. surface satisfying T e g = Tpo = 150 



FIG. 9. (Color online) The y/s dependance of the par¬ 
ticle production measure Aaniso/Mdeai — 1. Black squares 
and red dots denote the (2+l)D boost-invariant and (3+l)D 
anisotropic hydrodynamics, respectively. We compare them 
with the results obtained within (2+l)D boost-invariant vis¬ 
cous hydrodynamics (blue diamonds). 

MeV. In Fig. 9, blue diamonds and black squares present 
the calculation for boost-invariant versions of viscous 
and anisotropic hydrodynamic models, respectively. The 
(2+l)-dimensional viscous hydrodynamics results were 
generated using the code of Luzum and Romatschke 
[63]. We also show results for full (3+1 (-dimensional 
anisotropic code (red dots). We note that there is some 
quantitative uncertainty in the presented results due to 
the effective-temperature freeze-out prescription used. 
Another possibility for the freeze-out condition would be 
to use a constant value of the Knudsen number [64]. We 
have not considered this possibility in this work. 
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